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Abstract 

In this work, we analytically study the Schrodinger equation for the (non-pure) dipolar ion potential 
V(r) = qjr + DcosO/r'^, in the case of 2D systems using the separation of variables and the Mathieu 
equations for the angular part. We give the expressions of eigenenergies and eigenfunctions and study their 
dependence on the dipole moment D. Imposing the condition of reality on the energies En^m implies that 
the dipole moment must not exceed a maximum value otherwise the corresponding bound state disappears. 
We also find that the s states (m = 0) can no longer exist in the system as soon as the dipole term is present. 
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I. INTRODUCTION 


The solution of the Schrodinger equation is very important for the study of atoms and molecules. 
However the number of systems for which analytical solutions exist is very limited. This motivated 
the search for simplified models of these physical systems resulting in several studies on spherical 
potentials. Despite such simplifications only few potentials have been analytically solved. Besides, 
real physical systems, such as atoms and molecules, rarely exhibit simple spherical symmetry like 
that of the hydrogen atom. 


fl, 


and was then properly structured with the work of Hautot [^. Their results paved the way 


The study of non-central potentials began with the pioneering works of Makarov [l| and Hartman 


for realistic applications of non-central-potential problems such as ring-shaped organic molecules 
which include, for instance, cyclic polyenes and benzene [4|. Since then there has been a significant 
interest in the literature in studies of non-central and ring-shaped potentials (see for instance 


Refs. 0-7] and references therein). However, as demonstrated by Hautot only few of these 
potentials in fact have analytical solutions, and they have thus been dominantly considered with 
either numerical technics or approximation methods. 

Owing to the emergence of graphene 


0, 0 , two-dimensional systems have recently attracted a 
lot of interest in the domain of Material Sciences. For this reason, this paper is devoted to the 
analytical treatment of the dipole potential added to the Coulomb potential V(r) = a/r+ h cos 9/r"^ 
for 2D systems. This potential, which we call non-pure dipole, is applicable to the case of an ionised 
dipolar molecule such as water. It has been considered in the case of three-dimensional systems 


using perturbation method [ll|. Furthermore, the pure dipole potential bcosO/r"^ in 3D has been 
widely studied both in nuclear and molecular physics [l^, ^ (see also references in [i0b where it 
has been shown that the dipole moment must exceed a critical value in order to have bound states. 
On the other hand, in the case of molecular dipole or pure dipole potentials, the moment must be 
below a critical value in order for bound states to exist [l^. For 2D systems, the pure dipole was 


by AlHaidari using the approximation method of tri-diagonal matrices {lfl| and by Moumni et al 


studied in Refs. 


14l4lfi|| and it was found in Ref jl^ that the critical dipole moment for bound 


states to exist is zero. 

This paper is organised as follows. In section m we start by showing that this potential is 
the first non-central approximation when considering a non-zero charge distribution which is not 
spherically symmetric or when considering non-spherical ions. This potential is applicable for most 
physical systems excluding neutral molecules. We then solve the Schrodinger equation for this 
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potential, in section imi using the separation of variables. Finally we conclude by discussing our 
results in section m 


II. NON-PURE DIPOLE POTENTIAL 


Our aim is to solve the Schrodinger equation for a system consisting of a point charge q under the 
effect of an extended charge Q = qj (a cluster of point charges qj). The latter is characterized by 
a non-zero total charge and a non-spherically-symmetric distribution. One can take as an example 
of this system a polar ion and a point charge. The potential produced by the charge distribution 
at the position of the test charge q is written as follows 


V(r) = Tl^-- 

^ dvreo rj 


( 1 ) 


Here qj is the charge of the component of Q, and fj = AjM = OM — OAj = r — aj is the 
position of the point charge q relative to this component. We defined M as the position of the 
charge q (also denoted by the vector r) and Aj as that of qj (defined by the vector Sj) relative to 
the origin O, which we choose to coincide with the center of the charge Q. Thus we write 


F(r) = ^ 


1 


Qj 


dvreo \r — a A ^ dTreo 

I J ' -j 
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1 


Qj 


{r-djf 


- 1/2 




dvreo r 


1-2 


r ■ ai 


- 1/2 
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( 2 ) 


We assume that the dimensions of the extended charge Q are small compared to those of the 
whole system constituted by Q and the point charge q, such that we write \aj\ <C |F|, and thus we 
have 

r2\ 


r ■ a,- 

1 -2—J- 


dj\ ^ r ■ dj I CLj 

T 1 ~ 1 H——h (D I 

rp Z j p Z \ p Z 


(3) 


Taking into account these considerations and keeping only the terms up to order aj jr in the above 
expansion, one can easily write the potential as a multipolar expansion 

1 I Qi qjCij cos 9j 


Vir) = 


dvreo 


V — tV 


(4) 


3 3 

The first term in this expression is the Coulomb interaction between the total charge Q and the 
point charge g, or the monopole part, while the second is the effect of the geometry of the non- 
spherically-symmetric charge distribution Q, which represents a dipole part. We can simplify 


3 


















the dipole term by mapping the charge distribution Q into a single dipole whose poles are the two 
centers of all positive and negative charges contained within it. This leads to a dipole-ion potential, 


i.e. a non-pure dipole whose total charge is non-zero 


^ (Q QdcosO 
dvreo V’’ ^ 


IQ ID cos ( 
-:2L _|- 

dvreo r dvreo 


where d is a characteristic distance of the dipole, 6 is the angle between the position vector of q 
relative to the center of the dipole and the axis of this dipole, and D = Qd is the dipole moment. 


III. 2D SCHRODINGER EQUATION FOR A NON-PURE DIPOLE POTENTIAL 


The stationary Schrodinger equation reads 


-A + U ^ = E^. (6) 

2m \ 

Since we are working with a 2D system and referring to the shape of our potential, we use the 
polar coordinates 0 < r < oo and 0 < 0 < 27r, and thus we write the wave equation as follows 

^cosd \1 

2m \(9r2 ^ ^2 qq 2 j 47rgQ r ) 

This expression may then easily be written in the separate form 

f 8“^ Id 2mqQ 1\ 1/5^ 2mqD \1 2mE 

We use the separation of radial and polar variables to write the solution as ^{r, 9) = r~^^‘^R{r)Q{9), 
and to split the equation into two parts, angular and radial ones 

In this work, we use the same considerations as those for molecular systems; a positive extended 
charge and a negative point charge, which are equal in magnitude {q = Q)- Hence 


-^0 + \/2Dcos0^ 0(d) = 0, 


j^+[Ee + -]^ + - + E R{r) - 0 . 


For simplicity we choose q = — |e|, where e is the electric charge (i.e., the point charge is an 
electron) and, to reduce the length of the expressions, we employ the atomic Rydberg system of 
units where 2me = h = e^/2 = dvreo = 1- 
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In order to extract the energy eigenvalues of the system as well as the eigenfunctions 
we first solve the angular equation to find the eigenvalues Eg, which may then be used to solve the 
radial part. 


A. Solution of angular equation 


The angular equation can easily be 
substitutions 6 = 2z, a = —AEg and p = 


cast in the Mathieu equation form 

-2V2D 



, by making the 


d‘^e{z) 

dz'^ 


+ (u 


2pcos22;) 0(z) = 0 . 


( 11 ) 


Notice that since the period of 6 is 27r then that of 2 is vr, and hence the solutions to Eq. m are 


just the cosine-elliptic ce 2 m and the sine-elliptic se 2 m +2 solutions of the Mathieu equation, where m 


is a natural number 


18l |. According to Floquet’s theorem [l9| (or Bloch’s theorem j^), for a given 


value of the dimensionless parameter p, the solutions are periodic only for specific values of the 
parameter a, which are known as the characteristic values. We denote these values corresponding 
to cosine solutions by a and those related to sine functions by b. 


Looking at equations ([7]) to (fTO]) one can see that, when the dipole effects disappear, i.e., L) —)• 0, 
our system has as a limiting behaviour the Coulomb system. We therefore require that our solutions 
must tend to the Coulomb solutions in this limit. This means that we must only keep the cosine 
solutions and discard the sine ones sine they are not valid solutions when T> —)• 0 for m = 0, where 
m is the orbital quantum number, as we shall show later. 


In general there is no analytical expression for these solutions and they are usually given either 
numerically or graphically (see Fig. [T|). However from these values one can extract the eigen- 
solutions of the angular equation which can then be used to obtain the radial ones. 

Referring to the expression of the dipole potential ([5]), we see that the dipole term is just a cor¬ 
rection to the mono-polar term since it is just the latter multiplied by the factor a/r. Furthermore, 
since it depends on the dipole moment, then p can be treated as a small parameter and one can 
write p <C 1. In this limit, the characteristic values have asymptotic analytical expressions which 
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a2m(p) 



FIG. 1. a 2 m(p) for m = 0,1, 2, 3,4 from bottom to top (in a.u.). 


we write as a series in p [1 
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(12a) 

(12b) 

(12c) 

(12d) 


For values of the index of a 2 m greater than 6 (i.e., for m > 3), one can use the same expression for 
both the a’s and the 6’s, which include up to 0{p^) terms 


02m = 4m^ + 


1 

2(4m2 


2 20m^ + 7 4 

1) ^ "^32 (4m2 - 1)® (4m2 - 4) ^ 

+_36m^ + 232m^ + 29_ , ^ ^ ^ « 

64 (4m2 - 1)^ (4m2 - 4) (4m2 - 9) 


(13) 


Using the expressions of the characteristic values and the relations a = —AEq and p = —2y/2D, we 
get the eigenvalues Eg as a function of the electric moment of the system D 
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For m > 3 we have 


® (4m2 — 1) 


- 


20m2 + 7 


2 (4m2 — 1)'* (4m2 — 4) 


D^+ 


32 (36m^ + 232m2 + 29) 


■D^ + 0{D^). (15) 


(4m2 — 1)® (4m2 — 4) (4m2 — 9) 

We see that in the limit p —>■ 0 (or D —)• 0), the Mathieu equation has as solutions cosm9/2, with 
a = 4m^. Thus the characteristic values can in all cases be written as 


a2m = 4m^ + Pm{p) ■ (16) 

Similarly for the angular eigenvalues we have 

4'™) = + Pm{D) , (17) 

where Pm{p) and Pm{D) are polynomials which are written in terms of even powers of p and D 
starting from 2. 

Next we use the expression of 4^”^^ in order to solve the radial equation. 


B. Solution of the radial equation 


We rewrite the radial equation (jlObh in terms of 4^^^ as 


+E + -+(Ef^^+ ^ 


R{r) = 0. 


(18) 


In order to simplify Eq. (jlSp into a class of known differential equations we use the following 
ansatz: 


R (r) = r^e f (r) , 


(19) 


where A is a constant to be determined. By means of the substitution (jl9p . the differential equation 
for / (r) becomes 


+ 2 (A - /3r) — + 2 (1 - A/3) 


/ = 0 , 


where /3 and A satisfy the following relations 


P^ = -E, and A(A-1)+ 4^”"^ + ^ = 0. 


Solving for A in the latter equation yields two solutions 


A=i±yr^ 


( 20 ) 


( 21 ) 


( 22 ) 
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However since we require R (r) to be a nonsingular function at r = 0, then the accepted value of 
A is 


A = i + 


(23) 


Now taking z = 2/3r and substituting into Eq. (l20t) . the latter reduces to a differential equation 
of the confluent hyper geometric type 


d^ , . d ( 1 


f{z) = 0. 


(24) 


The solution of this differential equation which is regular at the origin 2 : = 0 is given in terms of 
confluent hypergeometric functions as 


f{z) = NiF,{X-/3-\2X,z), 


(25) 


with N a normalization constant to be determined later. 

In terms of the variables r and 9, we can now write the general form of the wave function -0 as 
follows 


(r, 9) = Nr^-^ 0 (9) iFi (A - 2A, 2/3r) . 


(26) 


The functions iFi (—n^, 2A, 2/Sr) may be written as Laguerre polynomials of degree as follows 

-1) (2/3r) = 1^1 i-nr, 2A, 2/3r) . (27) 

To determine the normalization constant N we substitute the wavefunction (|26p into the nor¬ 
malization condition f {i/j {r,9)\‘^ rdr d9 = 1, where we recall that 0 (9) is the Mathieu solution 


ce2m{9l2) which is normalized by definition [1^. Using the identity 


m 


f 


e-^ 




n 


(28) 


we obtain 


N = 


2A/JA+I 

(2A-1)! 


(n-k2A- 1)! 


(29) 


n! (n + A) 

Using the condition of the convergence of the solutions at infinity, and from the asymptotic 
behavior of the confluent series iFi of Eq. (1261) . that for r —>■ 00 we have lEi —)• 0 which leads to 
V’ —?■ 0 at infinity, we obtain the following general quantum condition 


A — /3 ^ = —Hr 


— 0 , 1 , 2 ,... , 


(30) 














and we get the discrete energy levels from the condition (13011 as follows 


Enr,m = -[nr + \J-E);"^’ + - 


d2m) 1 
2 


-2 


(31) 


We can relate these energies with the Coulomb energy by using the relationship (jl7p 


Enr,m — 




-2 


(32) 


and the limit Pm{D) —)• 0 gives us 


Enr,m = - inr + \m\- 


Pm{D) ^ 1 


-2 


2 2 


-2 


(33) 


which are the energies for the Coulomb potential with n = + |m|, where we have the condition 

m < n from the Coulomb system 2j, [23| . The final expression for the energy eigenvalues is 

-2 


En.m = - n - m + \/ -E, 


^(2m) 


1 

+ 2 


(34) 


The requirement that these energy eigenvalues -En,m be real necessitates that be negative, 

and thus from the relation a = —AEq the characteristic values must be positive. This gives 
us a condition on the parameter p and therefore on the dipole moment D from the definition 
p = —2y/2D. Consequently, we conclude that in order for bound states £'n,m to exist, it is necessary 
that the dipole moment D not exceed a critical value determined by the equation = 0. This 

critical value depends only on m and is thus denoted We show in table H] the values 44^™^ 

Rydberg atomic units) for some values of m. We note that when the dipole moment is non-zero 


m 

0 

1 

2 

3 

4 

5 

6 

7 

7^(^) 

-^crit 

0,000 

7,530 

24, 547 

51,285 

87,746 

133,930 

189,837 

255,468 


TABLE I. Critical value for some values of m. 


then the level £'o,o disappears. We shall elaborate on this point separately in the next section. 

From the relation (j3ip . we can plot all the graphs for the energies for all possible values of the 
dipole moment En,m{D), as shown in figures El El [H and [SJ 


IV. CONCLUSIONS 

In this work, we have analytically studied the Schrodinger equation for the potential V{r) = 
ajr + bcos9/r‘^ for 2D systems using the separation of variables and the Mathieu equations for 
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FIG. 2. En.i{D) for n = 1,2, 3,4, 5 from bottom to top (in a.u.). 
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FIG. 3. En^ 2 {D) for n = 2, 3,4, 5,6, 7 from bottom to top (in a.u.). 


the angular part. We have shown that this potential is applicable to a dipolar ion, and that it is 
the first-order approximation resulting from the effect of a non-spherical distribution of charges. 
We gave the expressions of eigenenergies and eigenfunctions and we studied their dependence on 
to the dipole moment D. As expected, our solutions tend to the Coulomb ones when the dipole 
moment vanishes. 

The plots of the energies show that they increase with the dipole moment up to a maximum value 
and then start decreasing. The behavior of these solutions is similar to that of the characteristic 
values of the Mathieu functions. Furthermore, the requirement that the energies F^n,m be real 
implies that the dipole moment must not exceed a maximum value, otherwise the corresponding 
bound states disappear. These critical values depend only on the magnetic quantum number m, 
and are thus labeled This result is similar to that found in 3D case by AlHaidari [l^. It is 

also in contrast to the case of a pure dipole (i.e., when no Coulomb term is present), where it is 
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FIG. 4. En^ 3 {D) for n = 3,4, 5,6, 7 from bottom to top (in a.u.). 





FIG. 5. En, 4 {D) for n = 4, 5,6, 7 from bottom to top (in a.u.). 


necessary that the dipole moment exceeds a minimum value in order for bound states to exist [Ij 
We also found that the s-states (m = 0) no longer exist for this system when the dipole term 
is present since = 0. Noting that the s-states exist for other non-central potentials (see for 
instance Refs. 7| and references therein), we deduce that the absence of these states in our case 
is not due to the difference between the symmetry of the s-states (which are central) and that of 
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the potential (which is not spherically symmetric). This phenomena requires further studies given 


the specihcity of the potential cos0/r^ 


M 


25l |. From a purely mathematical point of view, this is 


due to the fact that the characteristic value oq is negative for all possible values of the parameter 
p, while the other parameters am^o start from positive values and end in the negative domain. 
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